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ABSTRACT 

There is a growing consensus that in the present universe most baryons reside 
in galaxy clusters and groups in the form of highly ionized gas at temperatures 
of 10^ ~ 10^ K. The H-like and He-like ions of the heavy elements can produce 
absorption features - the so-called "X-ray Forest" - in the X-ray spectrum of a 
background quasar. We investigate the distribution of the X-ray absorption lines 
produced by this gas under three different cosmological models: the standard 
CDM with fio = 1, a flat model with Qq = 0.3 and an open model with with Qq = 
0.3. We give a semi-analytic calculation of the X-ray forest distribution based 
on the Press-Schechter formalism, following Perna & Loeb (1998). We choose 
three ions (O VIII, Si XIV and Fe XXV) and calculate the distribution functions, 
the number of absorbers along the line-of-sight (LOS) to a distant quasar vs. 
redshift and column density in a given ion. We find that significant differences in 
the evolution of the distribution functions among the three cosmological models. 
Using Monte Carlo simulations, we simulate the distribution of X-ray absorption 
lines for 10,000 random LOS. We find there are at least several O VIII lines 
with column density higher than lO^^cm^^. Finally we explore the possibility 
of detecting the X-ray forest with current and upcoming X-ray missions and we 
present an XMM RGS simulation of a representative quasar X-ray spectrum. 

Subject headings: galaxies: clusters: general — intergalactic medium — 
large-scale structure of universe — quasars: absorption lines — X-rays : general 
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1. Introduction 

There is an apparent deficit in the total density of baryons at moderate and low 
redshift, z^ which has come to be called the "missing baryon" problem (see, e.g., Fukugita, 
Hogan & Peebles 1998; Cen & Ostriker 1999). Observations of the hydrogen and helium 
absorption lines in the Lyman alpha forest give a baryon density at high redshift (z ~ 3) 
of 1^5 > 0.017/i^^ (Ranch et al. 1998), which is consistent with the value fit = 0.019/;,^^ 
derived from standard big bang nucleosynthesis (Buries & Tytler 1998, a Hubble constant 
of Hq = lOOh km s~^ Mpc~^ is used throughout the paper). However, in the local universe 
the baryon budget is far below this number. Summing over all the baryons inside stars, 
neutral atomic gas and molecular gas gives Qb ~ O.OOSh^^ . 

There is a growing consensus that these "missing baryons" reside in a hot (> 10^ 
K), ionized plasma associated with groups of galaxies. This is effectively attributing to 
the universe a mix of components similar to that observed in richer groups and galaxy 
clusters (Fukugita et al. 1998; Cen & Ostriker 1999). It is plausible that this medium can 
be enriched in heavy elements. These heavy elememts would not be fully ionized at the 
temperatures of interest, as is the cluster gas (Renzini 1997). 

Resonant absorption by a hot, enriched medium would introduce features in the 
X-ray spectrum of a distant quasar. The features would be narrow lines or broad troughs 
depending on the velocity structure of the absorber. The associated absorption edges would 
also be present but are generally much weaker (Markevitch 1999). Early work by Shapiro 
& Bahcall (1980) discusses the X-ray absorption spectrum introduced by a uniformly 
distributed, hot IGM with an admixture of metal atoms via X-ray "Gunn-Peterson" 
effect. Using the same method but with a photoionized model of the IGM, Aldcroft et al. 
(1994) constrain the density and temperature of the IGM with ROSAT PSPC spectra of 
z ~ 3 quasars. Basko, Romberg & Moskalenko (1981) discusses the detectability of X-ray 
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resonance absorption lines in quasar spectra produced by hot plasma in an intervening 
galaxy cluster. More detailed and accurate calculations of cluster absorption were performed 
by Gil'fanov, Syunyaev & Churazov (1987); Krolik & Raymond (1988); Sarazin (1989). 

Like the Lyman alpha forest system in the optical band, the spatial distribution 
of the galaxy groups and clusters can also produce an "X-ray forest" along the line of 
sight in the X-ray spectrum of a background quasar. The concept of an "X-ray forest" 
was first suggested by Hellsten, Gnedin & Miralda-Escude (1998): X-ray absorption lines 
are produced by hot intergalactic medium in the form of "filamentary and sheet-like 
structures connected to galaxy clusters and groups, as well as colder gas left out in voids." 
A similar concept was explored by Perna & Loeb (1998) using the expected spectrum of 
mass concentrations for a universe dominated by cold dark matter. The effect of differing 
cosmologies on the number and evolution of X-ray absorption lines is related to the effect 
on the number density of clusters (Bahcall & Fan 1998; Eke et al. 1998). This means 
that X-ray absorption line studies might eventually provide independent constraints on 
cosmological parameters. 

Here we build on and extend the approach of (Perna & Loeb 1998) to explore the X-ray 
forest for various cosmologies and to assess its detectability. Detecting it is not easy. Most 
of the absorption lines will have equivalent widths < 1 eV, and none of the previous X-ray 
missions {Einstein, ASCA, ROSAT) had sufficient sensitivity to detect such features in a 
quasar spectrum. However, current and future missions give order-of-magnitude advances in 
sensitivity for the X-ray forest. For example, Chandra grating spectrometers have resolving 
powers of 1,000 around 1 keV (ASC Science Center 1997) with sufficient collecting area 
to detect an absorption line from an ion column density of approximately lO^^cm^^ in a 
plausible quasar spectrum. Assuming a moderate metal abundance, this column density 
implies a hydrogen column density of a modest galaxy cluster (Canizares & Fang 1998). 
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XMM and Constellation- X achieve comparable or better energy resolutions with larger 
effective area. In this paper we use simulations to assess the ability of these missions to 
detect features in the X-ray forest. 

This paper is organized as follows: section II gives a semi-analytic calculation of the 
X-ray forest distribution function, based on Press-Schechter formalism. Section III is 
devoted to the numerical simulation of the distribution function. In section IV we discuss 
the detectability of the X-ray forest. Section V presents the conclusions and discussion. 

2. X-ray Forest Distribution Function 

2.1. Press-Schechter Formalism 

Galaxies, galaxy clusters and other large scale structures grow from the initial small 
scale density fluctuation via gravitional instability. The small scale fluctuation flrst grows 
linearly, until it reaches a critical density. Then it decouples from the Hubble expansion, 
starts collapsing and flnally condenses out as a virialized, gravitational bound halo. Given 
a random Gaussian distribution, the comoving number density of virialized halos can be 
described by the Press-Schechter function : 



1 



dn / 2 \ 2 p S^ 



dM„ir- \nj Ml- (r{z, M^i^) 



dlna 



''^^' 2a{z,M,„f^ ^^' 



d In M^. 

Here M^r is the mass of the virialized halo; p = SQqHq/StiG is the comoving mean 
density of the universe which is constant during matter domination; 6c denotes the linearly 
extroplated overdensity at which an object virializes; and (7{z,Mmr) is the rms density 
fluctuation inside halos containing a mean mass of Mmr- Press & Schechter (1974) flrst gave 
this function with a very simple and intuitive model : large scale virialized objects form 
from the nonlinear interaction of small scale objects through a self-similar condensation 
process. However, the original theory suffers the so-called "cloud-in-cloud" problem of 
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miscounting the underdense regions properly. Bond et al. (1991) and Lacey & Cole (1993) 
extended this model by counting the overdense regions one-by-one and gave the correct 
normalization of the mass function. Although no direct observational evidence shows 
that the Press- Schechter function is the right way to describe the cluster abundance, this 
function fits N-body simulations extremely well (Lacey & Cole 1994; Eke, Cole & Frenk 
1996; Tormen 1998; Frenk et al. 1999). 

Formally one would expect 5c to depend on the cosmological model and the geometry 
of the collapsing object. Since most of the rich clusters are fairly round it would be a good 
assumption that the collapse is close to spherical. In a flat universe spherical collapse gives 
5c = 1.686 (Lacey & Cole 1993). However, the value of 5c changes by < 5% as one goes 
from an Einstein-de Sitter universe to an ^2 = 0.1 universe (Eke, Cole & Frenk 1996). So 
we adopt the value of 1.686 throughout the paper. 

The rms mass variance at redshift z can be expressed by the present rms mass variance 
and the linear growth factor D{z) through 

a{z,M) = a{z = 0,M)^^^ (2) 

From observations of cluster density in the local universe, several papers give the 
normalization of the power spectrum at Sh^^Mpc scale (Eke, Cole & Frenk 1996; Viana 
& Liddle 1996; Pen 1998; Viana & Liddle 1999). Here we adopt the value from Eke, Cole 
& Frenk (1996). The present rms mass fluctuation o"(z = 0,M) can be calculated by the 
normalized power spectrum. A functional flt to D{z) is given by Lahav et al. (1991) 

To apply the Press-Schechter formalism we need to determine the virial mass of clusters 
(Myjr) precisely. However, measuring the virial mass of clusters is difficult, especially 
at high redshift. Observationally the X-ray temperature (Tx) of galaxy clusters can be 
measured more accurately. Numerical simulations and observations of X-ray clusters show 
that there exists a very tight relationship between Tx and M^ir (Hjorth, Oukbir & Kampen 
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1998; Bryan & Norman 1998). To simply we assume a singular isothermal sphere model of 
the cluster mass density distribution. The mass-temperature relationship is given by Eke, 
Cole & Frenk (1996) 

Here A^ is the ratio of the mean cluster density to the critical density at that redshift, (3 
is the ratio of the specific kinetic energy to thermal energy, VLq and VL{z) are the cosmology 
density parameter at present and redshift of z, respectively. 

Recently various authors show that equation (3) is accurately obeyed in N-body 
hydrodynamic simulations with value of /9 ~ 1 (Navarro, Frenk & White 1995; Evrard, 
Metzler & Navarro 1996; Bryan & Norman 1998). Although all of them suggest a slightly 
higher /5-value, it might be due to the incomplete thermalization of the intracluster gas, 
or the gas density dropping faster than r~^ around the virial radius in the numerical 
simulations. Here we use (3 = 1, which means that the specific galaxy kinetic energy equals 
the specific gas thermal energy within the virial radius. A Recent analysis on Abell 401 
shows the cluster mass given by the best-fit model is approximately by a factor of 1.7 lower 
than the value predicted by equation (3) (Nevalainen, Markevitch & Forman 1999). The 
difference is attributed to the fact that the simulated clusters have steeper gas density and 
shallower temperature profiles than observed. 



2.2. Gas Column Density Profile within Galaxy Clusters 

Assuming a "/3 model" of the cluster gas density distribution (Sarazin 1988), the 
column density of gas particles at an impact distance of b is (Perna & Loeb 1998) 



N (b) = faaskTx 



b-' 



(4) 



Here b is the projected distance from the center of the galaxy cluster; Vc is the core radius 
of the galaxy cluster and we select a constant value of 250 kpc throughout the paper; fgas 
is the baryonic gas fraction; fi = 0.59 is the mean atomic weight and k is the Boltzmann 
constant. 

We are interested in the metal ion column density which would produced absorption 
features in the X-ray spectrum of a background quasar. The ionization sources of 
intracluster gas can be either photoionization from the X-ray background radiation or 
collisional ionization. However the X-ray background is in general too weak to be the main 
source of ionization so we only consider collisional ionization here. Since generally the 
collisional time scale is much shorter than Hubble time, the gas is in collisional ionization 
equilibrium. If we denote T(X*) = N{X^)/N{X) as the fraction of ion X*, T would be only 
a function of temperature, i.e., T = T(T) under collisional equilibrium (Sarazin & Bahcall 
1977). Assuming a uniform metallicity Z{X) = N{X)/N{H) the ion X* column density 
distribution is 



N{X') =0.46 Z{X)T(T)- 



fgaskTx 



2G'rc(yump)2 
Here 0.46 is the fraction of hydrogen atoms by number. 



1 



h-' 



(5) 



The baryonic gas fraction within galaxy clusters has recently received attention (White 
et al. 1993; White & Fabian 1995). Almost all the observations give large values of baryon 
density Vti, than expected from the theory of the big bang nucleosynthesis if VLq = 1, which is 
taken as an indication of a low density universe. Both local and high redshift observations 
(Rines et al. 1999; Ettori & Fabian 1999) suggest fgas scatters between 0.1 and 0.3. To 
simplify, we use fgas = 0.2 throughout the paper. 
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2.3. Distribution Function 

In analogy to the Lyman-a forest system, we define d'^P/dN^dz as the number of 

absorption systems along the line of sight with a column density between A^* and A^* + dN^ 

per unit redshift. Here A^* = N{X'^) is the column density of ion X*. If we define S as the 

cross section of a galaxy group or cluster, the distribution function is then 

d'^P f ,_ dn di: de 

dTx---r-r— (6) 



dN^dz Jt^ dTx dN^ dz 
Here dn/dTx is the cluster number density distribution at different cluster temperatures, 
which can be obtained by the Press-Schechter function and the cluster mass-to-temperature 
relationship, i is the path length. We refer to Perna & Loeb (1998) for detailed calculations. 

3. Computed Distribution 

3.1. the X-ray Forest in Different Cosmological Models 

Given different cosmological models, the X-ray forest distribution d'^P/dN^dz is 
determined through equation (6) via the Press-Schechter function (equation (1)). Here we 
compare the X-ray forest distribution under three cosmological models: a standard cold 
dark matter model (SCDM); a low-density open CDM model (OCDM) and a low-density 
A-dominated CDM model. Table 1 shows all the model parameters. 

EDITOR: PLACE TABLE 1 HERE. 

Observationally resonance absorption lines will be produced by transitions with the 
ions of the most abundant elements with the largest oscillator strength. For specificity we 
choose the ions O VIII, Si XIV and Fe XXV. Table 2 gives all the parameters. Mushotzky 
et al. (1996) measured the metal abundances in four rich clusters with ASCA. We adopt 
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their metallicity because we concentrate only on high column densities produced by rich 
clusters. The redshift evolution of cluster metallicities can also play an important role here, 
for simplicity we treat them as constants. We adopt the solar abundances from Anders & 
Grevesse (1989). The atomic data are from Verner, Verner & Ferland (1996). Figure 1 
shows the ionization fractions (Mazzotta et al. 1998), in which O VIII, Si XIV and Fe XXV 
show three different peak temperatures around 2 x 10^, 10^ and 5 x 10^, respectively. These 
temperatures roughly correspond to the galaxy cluster mass around 10^'^Mq, IO^^Mq and 
IO^^Mq, according to the mass-temperature relationship (3). 

EDITOR: PLACE TABLE 2 HERE. 

EDITOR: PLACE FIGURE 1 HERE. 

To compute the Press-Schechter function we use the codes based on Eisenstein & Hu 
(1998) for the calculation of power spectrum and mass variance o"(M). However, their codes 
are normalized by COBE data, which is on a scale larger than the size of galaxy clusters, 
so we renormalize the power spectrum based on the observations of the cluster abundance 
in the local universe. 

Figure 2 through figure 4 give the distribution of three elements: Fe XXV, Si XIV and 
O VIII under three cosmological models at different redshifts. The most striking feature is 
that in all three ions, the standard cold dark model shows a more rapid evolution than the 
two other low-density models. This trend in general agrees with the evolutionary scenarios 
predicted by different cosmological models. In the SCDM model, the density perturbation 
grows as (l + z)^^, most clusters are formed recently and the cluster number density declines 
very quickly as we move to high redshift. On the other hand, low-density models predict 
slower evolution. Figure 5 shows the comoving cluster abundances as a function of X-ray 
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temperature at different redshift predicted by the Press-Schechter distribution. Although 
LCDM also shows some decline of cluster abundances at high temperature, SCDM show the 
strongest evolution of number density at temperatures > 5 x 10^ K. Among all three ions, 
the distribution of Fe XXV shows the biggest decline (over three decades between z = 
and z = 3). The reason is that Fe XXV has the highest peak temperature (~ 5 x 10^ K) in 
figure 1, which corresponds to the biggest decline in the cluster abundances (figure 5). 

EDITOR: PLACE FIGURE 2 HERE. 
EDITOR: PLACE FIGURE 3 HERE. 
EDITOR: PLACE FIGURE 4 HERE. 

EDITOR: PLACE FIGURE 5 HERE. 

The second feature is that all three distributions show a rapid cutoff between lO^^cm^^ 
and lO^^cm^^. This is mainly because the maximum ion column density can not exceed the 
central column density of the galaxy clusters. The central ion column density is given by 
equation (5). For instance, the peak temperature of O VIII is around 2.4 x lO^K (figure 1), 
which roughly corresponds to a central O VIII column density of 6 x lO^^cm"^. This 
explains why figure 4 shows a sharp cutoff around this column density. 

The third feature is that for all three ions, the SCDM model predicts more absorbers 
along the line-of-sight than OCDM and LCDM. This is because the spatial density of 
clusters is very sensitive to the present matter density - Qq (equation (1)). Qo represents 
the overall amplitude of the density fiuctuation and changing of Qq will dramatically change 
the spatial density of virialized objects, the clusters of galaxies. This can be a method to 
determine Qq in the future from observations of the X-ray forest. 
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3.2. Monte Carlo Simulations 

To obtain statistics on tlie distribution of X-ray absorption lines, we carry out a series 
of Monte Carlo simulations. Assuming both absorber column density and redshift are 
independent random variables (M0ller & Jakobsen 1990), we define the Probability Density 
Functions (PDF) for each variable. 



min 

«<^'' - aL WW/' 

Here A is the total absorber number, given by the integration of the distribution 
function over both redshift and ion column density. We set a range for each variable as 
A^* G [10^^, 10^^] cm"^ and z G [0,3]. Using these PDFs we obtain column densities and 
redshifts for 10,000 randomly selected lines-of-sight. Figures 6 and 7 show the average 
cumulative distribution of the absorption line numbers vs. column density and redshift, 
respectively. Although Fe XXV is too scarce to give any statistical information, both figures 
do reflect the three features of the X-ray forest we noted before. In flgure 7 we include only 
the absorption lines with column density higher than lO^^cm"^ because this is the lowest 
column density which is detectable by Constellation X (see the next section). In this flgure 
we see that, compared to the LCDM and OCDM models, SCDM presents a larger number 
of absorption lines. For instance, SCDM shows that on average ten O VIII absorption lines 
up to z = 1, where OCDM and LCDM give only four and three lines, respectively. 

EDITOR: PLACE FIGURE 6 HERE. 



EDITOR: PLACE FIGURE 7 HERE. 



-13- 

Assuming a velocity dispersion of 6 ~ 300 km sec~^, the line optical depth is obtained 
by t{i') = A^V(z/). Here 0"(z/) is the absorption cross section at frequency u (Spitzer 1978). 
The line-of-sight transmission is defined as D = e~^. Figures 8 to 10 show the transmission 
of Fe XXV, Si XIX and O VIII under the three cosmological models, to 2; = 3. 

EDITOR: PLACE FIGURE 8 HERE. 



EDITOR: PLACE FIGURE 9 HERE. 



EDITOR: PLACE FIGURE 10 HERE. 



4. Detectability 

Generally most X-ray absorption lines we discuss here are narrow, unresolved lines. 
For a weak absorption line, the equivalent width Weq is given by (Spitzer 1978) 



E 



/;{l-e.p(-.„e-)} 



Where E is the line energy, b and c are velocity dispersion and light speed. Here tq is the 
optical depth at line center, tq = NsX/n^^'^b, where N is ion column density, A is wavelength 
and s is the Einstein absorption coefficient. 

The detectability of the X-ray forest is limited by the spectral resolving power and 
effective area of the spectrometer. A weak, unresolved resonance line has Weq <^ AE, 
where AE is the bin width, determined by the instrument resolving power R (i? = E/ AE). 
Suppose we use an instrument with resolving power R to observe an resonance line around 
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energy E (keV). The source spectrum has a continuum intensity of Fx (in units of photons 
cm^^ s^^ keV^^) around E. Given an observing time T, the minimum detectable equivalent 
width for an unresolved absorption line is 

Here (S/N) is the signal-to- noise ratio and ^e// is the effective area, and we assume 
negligible background. 

For illustration we choose a typical spectrum of an X-ray bright quasar with photon 
index F = 2.5 and flux 1.0 x 10^^^ ergs cm~^sec~^ between 0.1 and 2.4 keV. The Galactic 
column density is 5.0 x 10^° cm~^. Given this representative spectrum we calculate 
the minimum detectable equivalent width and column density of some ion species for a 
particular observation time. In Table 3, we list three instruments: Chandra LETG/HETG, 
XMM RGS and Constellation- X Calorimeter/Gratings. Assuming the absorption ions 
are located at z = 0.5, the equivalent width is calculated based on a S/N of 3 and an 
integration time of 100 ksec. Comparing this table with figure 6 and 7, we find that O VIII 
ion is the best candidate for all three instruments. Using XSPEC 10, we simulate this 
"representative" spectrum plus O VIII absorption lines from one realization of the LCDM 
model (Figure 10) on XMM RGS. Several tens of quasars with z > 1 have such spectrum 
or are even brighter, so we put the redshift of this quasar aX z = 1. The LCDM simulations 
give three absorption lines with z < 1 and N > lO^^cm^^. Table 4 lists the simulated line 
properties. Then we fit the simulated spectrum with a model only containing Galactic 
absorption and a single power law. The x^ plot of figure 11 clearly shows three absorption 
lines. Notice here line 2 and line 3 are only separated by approximately 3.7 eV. 

EDITOR: PLACE TABLE 3 HERE. 
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EDITOR: PLACE TABLE 4 HERE. 

EDITOR: PLACE FIGURE 11 HERE. 



5. Summary 

In this paper we use a semi- analytic method to investigate the X-ray forest. Following 
previous work (Perna & Loeb 1998) we calculate the X-ray distribution function of O VIII, 
Si XIV and Fe XXV under the Standard cold dark matter model and two low-density 
models. Fe XXV shows the more rapid evolution of SCDM, compared to OCDM and 
LCDM; These trends are milder for Si XIV and O VIII. Using Monte-Carlo simulation, we 
investigate the average distribution of the X-ray forest. We find SCDM model presents 
more absorption lines than OCDM and LCDM, which eventually might yield a method of 
determining Qq. We also select a typical spectrum of an X-ray bright, distant quasar to 
explore the detectability of the X-ray forest. We find for all three telescopes, there are at 
least several O VIII absorption lines detectable by Chandra, XMM and Constellation- X. 
This result is consistent with Perna & Loeb (1998) and Hellsten et al. (1998). 

The X-ray forest distribution function depends on the Press-Schechter distribution of 
galaxy clusters, the cluster mass-temperature relationship and the gas distribution inside 
galaxy clusters (equation (6)). In the following we discuss several important factors which 
can affect the result. 

The key element of the Press-Schechter function is the mass variance, a{M). The 
present mass variance a{z = 0, M) is calculated from the present power spectrum filtered 
through a top-hat window, and normalized by the the mass variance at 8h~^ Mpc, (Tg. 
The estimation of erg involves fitting the Press-Schechter function with the observed spatial 
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number density of local clusters Henry & Arnaud (1991); White et al. (1993); Viana & 
Liddle (1996, 1999). In this evaluation a crucial relationship is the mass-temperature 
relationship of equation (3). Although there are both numerical and observational evidences 
for this relationship, it has a well-known problem of the "recent-formation approximation": 
the clusters observed today formed just before we observe them. To resolve this problem 
two different methods (Lacey & Cole 1993; Sasaki 1994) were proposed to substitute the 
Press-Schechter function, both of which gave nearly the same result (Viana & Liddle 1996). 
Based on the mergering-halo formalism of Lacey & Cole (1993), Voit & Donahue (1998) 
derived a new Mmr — Tx relation and claimed that equation (3) overestimated temperature 
evolution and so the numbers of high-z clusters. Their conclusion would affect the X-ray 
forest distribution function of Fe XXV discussed here, but not O VIII. The reason is that the 
exponential term of Press-Schechter function will only become crucial at high temperature 
(> lO^K), well above the temperature of the ionization peak of O VIII at ~ 2 x lO^K. 
However it would be important to investigate this effect because it would largely decrease 
the possibility of detecting the highly-ionized heavy metal absorptions lines, such as 
Fe XXV. 

Another important uncertainty which can affect the X-ray distribution function is the 
metal abundance inside clusters. Using ASCA^ Mushotzky et al. (1996) shows the mean 
abundances of O, Si and Fe of four galaxy clusters are 0.48, 0.65 and 0.32 respectively, 
which are close to the values used in this paper. However, there are a few factors which can 
affect the metal distribution. First, the metal abundances can depend on the mass of the 
cluster or group. Recent observations by Hwang et al. (1999) show at temperature above 1 
keV, the metal abundances are roughly 0.3 Mq, with little variance, while at temperature 
below 1 keV, the metal abundances drop very fast (Renzini 1997; Davis, Mulchaey & 
Mushotsky 1999, and references there in). If it is real, this effect can dramatically drop 
the possibility of observing Oxygen absorption lines because of its low peak temperature of 
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ionization. Another factor is the assumption in this paper of constant metal abundances 
upto redshift as high as z ~ 3, ahhough this is important only for the richest systems which 
have significant column density. No direct evidence shows a constant metal abundances 
beyond z ~ 0.3 (Mushotzky & Loewenstein 1997). 

Further progress in this subject relies on both numerical simulations and observations. 
Large-scale simulations on the X-ray clusters would provide us more accurate information 
on the X-ray forest distribution. On the observation side, with the launch of Chandra 
and XMM, we would expect a few O VII or O VIII absorption lines by observing low and 
moderate redshift quasars. In the future, Constellation-X will provide us superior spectrum 
of the X-ray forest, with the pioneer of probing the number and distribution of objects in 
the universe. 

We are grateful to Greg Bryan for many useful discussions and helpful suggestion on 
calculating the power spectrum. We thank Rasalba Perna for kindly discussions on her 
paper. The code for calculating power spectrum is provided by Daniel Eisenstein and 
Wayne Hu. We would also like to thank MIT/CXC team for help on using X-ray softwares. 
We also thank the suggestions from the referee. This work is supported in part by NASA 
contract NAS 8-38249. 
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Fig. 1. — Ionization fractions for O VIII (solid curve), Si XIV (dashed curve) and Fe XXV 
(dot-dashed curve). 

Fig. 2. — The distribution functions of number of absorbers per redshift and column density 
for Fe XXV. For each of the three cosmological models, the distribution function vs. column 
density is given at four different redshifts from top to bottom : z = O.0(solid line), z = 
1.0(dotted line), z = 2.0(dashed line), z = 3.0(dot-dashed line). There is little evolution in 
the f2o < 1 cosmological models, comparing to a more rapid evolution in the fig = 1 model. 

Fig. 3. — The distribution function for Si XIV. Line symbols are the same as the previous 
figure. The evolution of f^o = 1 model is still faster than low density models, but is milder 
than Fe XXV. 

Fig. 4. — The distribution function for O VIII. Line symbols are the same as the previous 
figure. 

Fig. 5. — Temperature distribution of galaxy clusters at different redshifts, predicted by the 
Press-Schechter function. For each of the three cosmological models, temperature distribu- 
tion are plotted from z = to z = 3. SCDM shows a rapid evolution at high temperature, 
compared to OCDM and LCDM. 

Fig. 6. — Cumulative column density distribution of average absorption line number upto 
z = 3. Each plot gives the absorption line numbers for one ion under the three cosmological 
models. O VIII shows more absorption lines. For instance, in the SCDM model, O VIII gives 
twenty-nine lines with column density over lO^^cm^^, compared to five for Si XIV and one 
for Fe XXV. 
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Fig. 7. — Cumulative redshift distribution of average absorption line number upto z = 3. 
Each plot shows the absorption line number with column density higher than lO^^cm^^ under 
the three cosmological models. In general SCDM shows more rapid evolution. 

Fig. 8. — Transmission of Fe XXV under the three cosmological models. The line energy is 
6.70 keV and spreads over redshift space from z = to z = 3. 

Fig. 9. — Transmission of Si XIV under the three cosmological models. The line energy is 
2.01 keV and spreads over redshift space from z = to z = 3. 

Fig. 10. — Transmission of O VIII under the three cosmological models. The line energy is 
0.65 keV and spreads over redshift space from z = to z = 3. 

Fig. 11. — The upper plot is the simulated spectrum, containing three O VIII absorption 
line at 440.11 eV, 464.82 eV and 468.48 eV, respectively. We fit the simulated spectrum by 
a simple powerlaw plus the Galactic absorption. The lower plot shows the x^, which clearly 
indicates three absorption features at the corresponding energy. The four big gaps in the 
upper plot are due to the instrumental gaps. 



-24- 



Table 1. The Cosmological Models 

Model fio ^R ^A h"" n^ dg 

SCDM 1.0 0.0 0.0 0.7 1.0 0.52 
OCDM 0.3 0.7 0.0 0.7 1.0 0.87 
LCDM 0.3 0.0 0.7 0.7 1.0 0.93 

'^Here we set Hq = 70 km s^^Mpc^^ 
"^We choose a scale-invariant power spec- 
trum with index n = 1 



Table 2. Ion Transition Parameters 



Ion Energy (keV) Oscillator Strength Abundance*^ Solar Abundance^ 



VIII 


0.65 


0.416 


0.5 


8.53 X 10-4 


Si XIV 


2.01 


0.416 


0.5 


3.58 X 10-5 


FeXXV 


6.70 


0.798 


0.3 


3.23 X 10-5 



'^Relative to solar abundance Zq. 
'^Relative to hydrogen number density. 
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Table 3. Resonant X-ray Absorption Lines at z = 0.5 

Ion Instrument Equivalent Width (eV) Column Density (cm" 

O VIII Chandra (LETG) 0.69 3.5 lO^^ 

XMM (RGS) 0.21 7.6 lO^^ 



Constellation X (Gratings) 0.04 1.3 lO^^ 

Si XIV Chandra (MEG) 1.34 4.7 10^^ 

XMM (RGS) 2.25 1.1 lO" 

ConsteJMioii X (Calorimeter) 0.09 3.110^5 

Fe XXV Chandra (REG) 24.7 1.7 lO^" 

XMM (RGS) N/A^ N/A^ 

Constellation X (Calorimeter) 0.52 9.2 lO^^ 

''The energy range of XMM RGS is 0.35 - 2.5keV 



Table 4. Selected O VIII Absorption Lines for XMM Simulation 



line number Column Density (cm ^) Redshift Energy (eV) Equivalent Width (eV) 



1 
2 
3 



4.28 10^6 


0.486 


440.11 


0.54 


2.80 10^6 


0.396 


468.48 


0.80 


3.27 10^6 


0.407 


464.82 


0.79 
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